The asymptotic relative entropy of entanglement 
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We present an analytical formula for the asymptotic relative entropy of entanglement w.r.t. PPT states for 
Werner states of arbitrary dimension. We then demonstrate its validity using methods from convex optimization. 
This is the first case in which the asymptotic value of a subadditive entanglement measure has been calculated. 
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The systematic investigation of quantum entanglement is a 
major goal of quantum information theory [|l]]. In the study of 
entanglement manipulation one considers the interconversion 
of different forms of entanglement by means of local quantum 
operations and classical communication (LQCC). For pure bi- 
partite states entanglement manipulation in both the finite and 
asymptotic limits is quite well understood. For pure states 
necessary and sufficient conditions for the local interconvert- 
ibility of entangled states are known. In the asymptotic limit 
of infinitely many copies of a pure state, a single number, the 
von Neumann entropy of a subsystem, appropriately quanti- 
fies the degree of entanglement |Q| . 

Much less is known about the entanglement of mixed states. 
One approach is to define entanglement measures, which are 
functions of a state that cannot increase under local opera- 
tions and provide constraints on possible local entanglement 
manipulation protocols. These measures prove to be useful 
mathematical and conceptual tools, and have interesting links 
with other areas such as the study of channel capacities [|^]. 
A number of such measures have been proposed, most no- 
tably the entanglement of formation [^, the distillable entan- 
glement and the relative entropy of entanglement [0|5]|. 
The distillable entanglement is defined as the asymptotic num- 
ber of pure maximally entangled states that can be obtained 
via LQCC from a supply of a given state. For mixed states the 
distillable entanglement is exceedingly difficult to compute as 
it is defined as an asymptotic quantity referring to infinitely 
many copies of a quantum state. Therefore, upper bounds on 
the distillable entanglement, in particular other entanglement 
measures, are of major practical use. One such entanglement 
measure is the relative entropy of entanglement, defined as 

En{<j)^n^inS[a\\p), (1) 

for states a, where T) is the set of states with positive partial 
transpose (PPT states), and S{a\\p) = tr[cr Ig cr — cr Ig p] (Ig 
signifies log2). This function essentially quantifies the distin- 
guishability of a from the set of PPT states. The set V can 
also be taken to be the set of separable states [^j]. However, 
the set of PPT states is much easier to characterise, and the 
resulting measure provides a tighter bound to the distillable 
entanglement, one that is actually attained on pure states and 



certain mixed states [^l . 

In general, efficient protocols for entanglement manipula- 
tion require an asymptotic number of copies of the initial 
state. Therefore, to address any question related to asymp- 
totic entanglement manipulation, one will instead have to con- 
sider asymptotic versions of the entanglement measures. For 
a given measure of entanglement E, the asymptotic version 
E°° is defined as the average entanglement per copy in the 
asymptotic limit [^0], 

E^{a) = lim '-. (2) 

For example, the asymptotic cost of creating a mixed state by 
LQCC from a supply of pure maximally entangled states is 
given by the asymptotic entanglement of formation [^. How- 
ever, such asymptotic entanglement measures are difficult to 
compute, and so far this task has not been accomplished ex- 
cept for the very rare occasions where the measure in question 
is known to be additive 
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FIG. 1. The asymptotic relative entropy of entanglement 
Er (a"(p)) of Werner states in ® as a function of p (the weight 
of the antisymmetric state) for several values of d. The dashed lines 
separate the two regimes of Eq. (^, for d = 10 and for d — 3. 

In this Letter we present a general formula for the asymp- 
totic relative entropy of entanglement for an important class 
of bi-partite states, namely the Werner states of arbitrary di- 
mension [^]. These states are the only states that are invari- 
ant under local unitaries of the form U ® U . They can be 
conveniently formulated in terms of the projectors ttq (tti) on 
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the symmetric (antisymmetric) subspaces of a Hilbert space 
7i = (8) C^. Denoting the permutation operator that inter- 
changes the subsystems as tt, the projectors can be expressed 
as ttq = (1 + 7r)/2 and tti = (1 — 7r)/2. A general Werner 
state is of the form (t{p) = pai + (1 — p)cto, p G [0, 1], where 
Co = 7''o/tr[7ro] and cti = 7ri/tr[7ri]. Appropriately exploiting 
the symmetry of this state is one of the key ingredients to the 
proof of the main statement of this Letter. We formulate this 
main result as a theorem. 

Theorem. For states in (g) C'' of the form a{p) — pai + 
(1 -p)cro withp e (1/2, 1], i?^ w.r.t. PPT states is given by 



I -Hip), 



^ - 2(i 



(3) 



lg^ + (l-p)lgf-i, p>^ 



where H{p) = —p\gp — (1 — p) lg(l — p). 



It is interesting to note from this formula that never 
exceeds the logarithm of the negativity [[l0|]. Curiously, 
E'^{a{p)) is a convex function of p, whereas numerics show 
that Eii{a{p)'^")/n is not convex for any finite number n > 
1. It is also intriguing that this formula is exactly the same 
as that from a different optimization problem considered by 
Rains [[Tltl, where he minimizes the function B{p, a) = 
S{p\\(t) + logja^l over all states a. Rains' quantity is also 
an upper bound to the distillable entanglement. However, it is 
not a convex function, which makes its minimization difficult 
due to the possibility of local minima. This applies especially 
in the asymptotic limit as the additivity properties of B{p, a) 
are unknown. 

Proof: The proof proceeds in two stages. First we provide 
an upper bound on Ej^ {a®'^)/n for any number of copies n, 
by presenting a trial optimal state. Then we derive a lower 
bound for Ej^{a'^")/n using convex optimization methods. 
We conclude by showing that the two bounds coincide in the 
asymptotic limit n oo. 

Upper bound: We consider a situation where we hold n 
copies of a Werner state a{p). Following the work of Ref. 
JT^], we will make heavy use of symmetry. In Refs. 
it has been shown that if a state is invariant under a certain 
symmetry group, then one can restrict the minimization in Eq. 
(|l|) to those PPT states that are also invariant under the same 
group. As a{p) is invariant under the group G of local uni- 
taries of the form U ®U, the state cr(p)®" is invariant under 
G®". This implies that the optimal PPT state for (t(p)®" can 
be chosen to be a convex sum of the 2" possible n-fold tensor 
products of ctq and ai. 



2"-l 
f=0 



(4) 



where x — (xq, X2^^^i)'^ forms a probability distribution, 
and fi e {0, 1}, i — 1, n, is the i-th digit in the binary 
representation of /. Furthermore, as the state cr®" is invari- 
ant under any permutation of the n copies, we can add the 



constraint that xi = Xm if the number of I's in the binary 
representations of / and m are the same. 

We will consider trial states of the form J^i^ii'^i'^i + 
WiCTo)*^" , which automatically satisfy this constraint. Such 
states correspond to: 



(5) 



where the Ui{vi) component is the weight of cri(CTo), Ui+Vi = 
1 and J2i Wi ~ I (and therefore J^j ^ !)■ 

We will need to know the eigenvalues of the partial trans- 
pose of rj in order to ensure that it is PPT. Ignoring degeneracy, 
there are only two eigenvalues of the partial transpose of the 
Werner state (j{p). It is easy to show that they are non-negative 
iff the following two component vector is non-negative: 
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1 )[l-p 



T 



1-p 



(6) 



Similarly, it can easily be shown that r/ will be a PPT state 
iff the y and x in the following equation are non-negative vec- 
tors: 



y 



-1 



1 



1 + 



f = r®"f (7) 



Subject to this condition we need to calculate the relative en- 
tropy between the n-copy Werner state (t(p)®" and this gen- 
eralized Werner state ry. This is given by: 

S{<j{pf^H/n 

n 

= -H{p) {lln)Y^C';y'\l~p)HgY,w.ur''vl (8) 

fc=0 i 

At this point, we notice that since the second term in Eq. 
(M) is the average of the function 



^(n - fc) = (1/n) Ig V w^u 



(9) 



over a binomially distributed variable k, the value can be sub- 
stituted by ^ {np) when we take the limit n ^ oo [jTsjl . 

Since E'^{a{p)) is the minimal value of S{<7{p)^"-\\r])/n 
over all possible PPT states rj, any such PPT state 77 will give 
us an upper bound for E^ {a{p)). In particular, for the vector 
X we propose a mixture of two n-fold Kronecker powers: 



E 



Ul. 



1 - fli 

a.; 



(10) 



where < Wi < 1, u>i = 1 — W2, < ai < 1/2 and 
1/2 < 02 < 1. Proper values for the parameters need to 
be selected, to ensure that the corresponding state r] will be 
PPT. Inspired by numerical results, we consider two separate 
intervals for p in Eq. (|); 1/2 < p < (d + 2)/(2d), and 
{d + 2)/{2d) <p<l. 
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For the first interval, set wi = and 02 = 1/2, so that 
Xfc = 2~" for all k. This state gives us an upper bound that 
equals Eji for one copy of a{p) 

n 

Now consider the second interval. We will set 



ai = 



( q! + 2)(1-p) 
d + 2-Ap 



02 = 



1 + d- (rf + 2)ai 
d + 2 - 4ai 



Wl = 1 — 1/72 



1 



1 + Z^' 



z={d + 2-4ai)/d. (12) 



We calculate T®"f : 



1 ^ 

SO that, with the values of Eq. ( p^ 



(13) 



Vk 



(14) 



where #(fc) is the number of I's in the binary representation 
of k. It is easy to check using Eq. ( [T^ that yk is always non- 
negative. 

As a consequence, by taking into account the discussion af- 
ter Eq. (^ and the notation introduced in Eq. (^, the upper 
bound for n ^ 00 reads 



l/n 



E^{a{p))<-H{p)^ lim Ig K] 



pn (1— p)n 



-ff(p)-lglim (%^) 
max(ti. t2) 



n \ 1/" 



- Ig 



max(qi,92) 



(15) 



with 



= d{d - 2)P(d + 2)i-P(l - 
^2 = (d - 2 + ^2(1 _ p))P{d'^p -d- 2)1- 
(ji = - 4 and 92 = ^(rf + 2 - 4p). 



(16) 



It is easy to check that for p > {d + 2)/ {2d) (the second in- 
terval), both ti > t2 and gi > 92- So we obtain: 



E^{a{p)) < ~H{p) - Ig ^ Ig ^Li +plg 1±|^ 



tl7) 



In the other regime, p < (d + 2)/ {2d), the bound obtainable 
from Eq. ( p3[ ) is worse than Eq. (pi]). This ends the proof of 
the upper bound. 

Proof of lower bound: We now proceed to find a lower 
bound on Eji. To do this, we will use the idea of Lagrange 



duality [|_4|1. To calculate E'^ we need to solve the optimiza- 
tion problem 

= mm{— iifp) z ■ \gx\ 

n X n 



with f > , -T«^"x < 



2"-l 

E 

k=Q 



Xk = 1, 



(18) 



where = (p, 1 — p)®". This is a convex optimization prob- 
lem, so it is possible to consider the so-called dual problem. It 
is a basic result in convex optimization theory that the optimal 
(maximal) g-value of the dual problem is a lower bound on the 
optimal (minimal) value of the primal problem, which is just 
what we are looking for (see [|l4|] for a general description of 
duality in optimization). 

The dual problem can be obtained as follows. First form 
the Lagrangian by multiplying the constraints with Lagrange 
multipliers and adding them to the objective function: 
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C{x,X,v) = —H{p) — —z^ Igx 
n 



2"-l 



-\TT®^x^v{Y,Xk-\) 

k=0 



(19) 



The constraint x > is not included explicitly, as it just deter- 
mines the domain of the function Ig. Note that A must exhibit 
the same copy-permutation symmetry as x. Because the con- 
straint associated to A is an inequality, we have to introduce 
the constraint A > 0. The dual function is now given by 



g{\,iy) = inf C{x,\,v), 

x>0 

and the dual optimization problem is: 

max (7(A, ly) , A > 0, 



(20) 



(21) 



including any other constraints on the domain of g. 

For our problem, the dual function can be calculated explic- 
itly. The derivative of the Lagrangian w.rt. Xk is 



dC _ 1 Zk 
dxk n In 2 Xfc 



(22) 



The Lagrangian reaches an extremum (minimum) at 

Xfc = -^, (23) 
^fc rt In 2 

where = v ~ |^r^"A^ and we have exploited the sym- 
metry of T. The dual function is 

,V , l + ln(nln2) 1 , 
nln2 n ^-^ 

k=0 

where we have used that ^ = {p,l ~ p)®", which implies 
that J2k^^ ■^k Ig Zk = —nH{p) . As stated before, A must be 
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non-negative, and inspecting the domain of g{\, v) yields an 
additional constraint that the be non-negative. 

Now, any acceptable assignment of values to A and v gives 
a lower bound to E'^{a{p)). Again we consider the two p- 
intervals of Eq. (|]). Forp < (d + 2)/ {2d) we propose 



= l/(nln2) 



d+1 



(8n 



d+l-2dp 
d+1 



(25) 



After a short calculation we obtain /i^. = i/^2"p" #('=)(! — 
p)'^{k) . This gives a feasible point of the dual problem, be- 
cause both A and as given here are non-negative. The value 
of g using these assignments is .g(A^, v^) = 1 — H{p). 

For the second interval, p>p' = {d + 2)/ {2d), we replace 
p by p' in Eq. (ESb, giving 



- , A^ = 




-1 
d+1 



(26) 



d+2 ■ 



We now obtain ^(A^, „^) = \g + (1 ^ p)lg !t. 

As the two lower bounds g{X^, i^^) and g{X^, v^) coincide 
with the two upper bounds (Eqs. (jllj) and (|l7|)), the proof of 
the Theorem is now complete. r-i 

The remarkable behavior of E^{a{p)) is shown in Fig. ll^ 
for several values of d. The nonlinear behavior for small val- 
ues of p goes over into a linear dependence on p above the 
threshold value p' = {d + 2)/ 2d. An immediate consequence 
of the result is that there are no inseparable Werner states with 
zero entanglement cost — a similar conclusion could not be 
drawn from Rains' bound as it is not an asymptotic quantity. 
It is astonishing that, as long as p < p' , is invariant un- 
der the strongly irreversible operation of twirling, mapping 
Werner states on to Werner states on ® 

Interestingly, the dependence of E]f{a{p)) on p is quite 
similar to the conjectured behavior of the entanglement of for- 
mation for a single copy of an isotropic state [p^: there, one 
can also distinguish between two regimes, and for larger val- 
ues of the weight F of the maximally entangled state in the 
isotropic state the dependence of the entanglement of forma- 
tion is conjectured to be linearly dependent on F. 

In this Letter we have concentrated on the important class 
of Werner states. With similar methods, one can also inves- 
tigate other classes with high symmetry. It is hoped that this 
work can significantly contribute to the quest for a better un- 
derstanding of the asymptotic regime of entanglement. 
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